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We present the theory of a Josephson parametric amplifier employing two pump sources. Our 
calculations are based on Input-Output Theory, and can easily be generalized to any coupled sys- 
tem involving parametric interactions. We analyze the operation of the device, taking into account 
the feedback introduced by the reaction of the signal and noise on the pump power, and in this 
framework, compute the response functions of interest - signal and idler gains, internal gain of the 
amplifier, and self-oscillation signal amplitude. To account for this back-action between signal and 
pump, we adopt a mean-field approach and self-consistently explore the boundary between amplifi- 
cation and self-oscillation. The coincidence of bifurcation and self-oscillation thresholds reveals that 
the origin of coherent emission of the amplifier lies in the multi-wave mixing of the noise components. 
Incorporation of the back-action leads the system to exhibit hysteresis, dependent on parameters 
like temperature and detuning from resonance. Our analysis also shows that the resonance condition 
itself changes in the presence of back-action and this can be understood in terms of the change in 
plasma frequency of the junction. The potential of the double pump amplifier for quantum-limited 
measurements and as a squeezer is also discussed. 

PACS numbers: 42.60.Da, 42.65.Lm, 52.35.Mw, 84.30.Le, 84.40. Dc, 85.25.Cp 



I. INTRODUCTION 

Recent advances in nanoelectromechanical systems (NEMS)^212I anc j optical communication technology^! have 
renewed interest in the field of quantum-limited amplification and encouraged scientists working in these areas to 
design systems which can operate at the quantum limit. The seminal work of Caves^ set the minimum noise added 
by a phase-preserving amplifier at half a photon referred to the input channel. The same work showed, however, 
that phase-sensitive amplifiers are not submitted to this limitation and can amplify one of the quadratures noiselessly 
at the expense of de amplification of the conjugate quadrature. Parametric amplifiers where a strong pump wave at 
frequency uj p causes simultaneous generation of signal and idler photons at frequencies u>s and u>j respectively, are 
particularly promising candidates to reach the quantum limit. These may employ either a second-order nonlinear 
interaction with lo p = ujs + u>i or a th ird-order interaction with 2uj p = lus + luj . Various schemes based on quantum 
optical parametric amplifie nJ 7 * 8 * 9 * 10 * 1 -^ operating at sub-quantum noise limits have been tested successfully. 

Low noise amplifiers in the RF and microwave domains have seen a growing demand from applications in astro- 
particle physics^E^ and quantum information processing^ involving microwave investigation of solid state qubits. 
Parametric amplifiers based on Josephson tunnel junctions (JTJ) are a front-runner for such applications as JTJs 
are the only known non-dissipative and nonlinear circuit elements, operating at radio frequencies and arbitrarily low 
temperatures. Experimental efforts^ 16 * 17 * 18 * 1 ^ ^ have engineered the nonlinearity of JTJ and utilized it to perform 
challenging measurements like non-demolition readout of superconducting qubitsSH. 

A typical parametric amplification process is depicted in Fig. [I] The amplifier can be visualized as a black-box 
where an incoming signal (A™) is reflected into a larger signal due to energy transfer from the pump drive to the 
signal. Similarly the idler port, empty of incoming waves, gives rise to a spontaneously emitted signal. The interaction 
between the signal and pump waves can be understood as a coupling between the internal coordinates of the amplifier 
representing the signal/idler (A) and pump (Y) modes respectively. Under the stiff pump approximation, the internal 
coordinate Y is considered to be completely enslaved to the pump amplitude Ap. As the noise at the input grows, 
one eventually enters the soft pump regime where the pump coordinate Y has its dynamics partly determined by X. 
This interaction between the pump and the signal leads to a reciprocal effect of the signal on the pump - referred to as 
'back-action' in the following text. As a result of this back-action, the outgoing pump wave (A P ut ) suffers depletion, 
which is the counterpart of amplification of the signal and idler waves. We have analyzed the soft pump regime and 
the effect of two-way coupling between signal/idler and pump channels in a detailed manner for a parametric amplifier 
using a single tunnel junction by using a model employing two pump sources. Interestingly, the use of dual-pump 
scheme simplifies the dynamics of this nonlinear coupled system significantly, a reasonable trade-off at the cost of a 
second generator. The absence of this crucial simplification makes a similar soft pump analysis for conventional single 
pump devices unwieldy.™ 

The quantum optics community is well-acquainted with the use of dual pumps in parametric amplification^ ^ 3 ! 24 !. 
The Josephson dual- pump amplifier (DPA) , besides providing an efficient theoretical tool for studying various aspects 
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AMPLIFICATION DEPLETION 

FIG. 1: Scattering representation of a parametric amplifier with its signal, idler and pump ports as well as its internal coordinates 
X and Y. The bidirectional coupling between the internal coordinate at signal frequency (X) and that oscillating at the pump 
frequency (Y) is shown using the curved arrows. This coupling leads, in one direction (left arrow) to amplification of the 
outgoing wave signal while in the other direction (right arrow) to depletion of the outgoing pump wave. It may be noted that 
this process also leads to a deterministic signal of finite amplitude at the idler port as the pump also populates the outgoing 
idler wave along with the signal, even if there is no idler input. 




of the dynamics of the amplification process, also offers practical advantages like greater bandwidth 2 -, tunable band 
center and additional degrees of design freedom. In the following sections, we provide a theoretical study of this novel 
Joscphson "paramp" and also explore the boundary between amplification and delicate effects like self-oscillation. 

In a regular degenerate paramp the signal is located within a cavity lincwidth of the pump frequency which is not 
favorable for the observation of fragile effects like self-oscillation. With a double pump scheme, the pump frequencies 
lie towards outer edges of the band of amplification, leading to better separation of the pump and the signal frequencies 
(Fig. |2j. This provides a zero background for detection of a small signal like self-oscillation which would have been 
otherwise challenging due to the large background of the pump. The use of dual pumping further helps in the 
realization of a symmetric phase boundary for bifurcation. This should be contrasted with a single pump device 
like the Josephson Bifurcation Amplifier (JBA)^, where the phase boundary is canted due to a first-order transition 
between dynamical states of the junction near threshold (cf. Appendix |A|. As a consequence, unlike the JBA, there 
is no inherent frustration between the thresholds for self-oscillation and parametric gain in the double-pump amplifier 
discussed here. It is this symmetry, unique to the DPA, that proves extremely valuable for the simplification of 
equations governing the dynamics of the system near bifurcation. We deal with this issue in greater detail in section 

This article is organized as follows: we first derive the basic equations governing the dynamics of the DPA in 
section |TT] showing the presence of various couplings existing in the system. Then in section |III[ we calculate the 
various response functions of the system for a stiff pump. In particular, we study the response of the system to noise 



and demonstrate the proficiency of the amplifier for achieving the quantum limited noise temperature in III A This 
is followed by a derivation of the steady state response of the system (section |lll B I where we show that the system 
exhibits self-oscillation beyond the bifurcation threshold. Also, we describe the use of the DPA as an efficient squeezer 
in section III C In the second part of the paper, we present a complete solution of the problem, accounting for the 
depletion of the pump due to back-action. We rederive the signal and self-oscillation amplitudes in the presence of 
back-action in section |TV} using a mean-field approach. In section [VJ we arrive at explicit expressions quantifying the 
magnitude of back-action in terms of system variables. Using the derived value of back-action we plot the corrected 
values of various response functions of the system. Finally we discuss the implications of our results and offer future 
perspectives in section IVI] 



II. BASIC MODEL OF THE DPA 



A circuit realization of the DPA is depicted in Fig. [3j It consists of a nonlinear LC oscillator pumped by two 
RF currents denoted by /^^(t) and Inp 2 (t). The nonlinear inductance of the circuit is provided by a Josephson 

tunnel junction which can be modeled as a linear inductance Lj — A- in series with a current-dependent inductance 

,2 

SLj = \Ljyi. Here we neglect higher order nonlinear terms since, for simplicity, the amplifier is operated in the 

jo 

weakly nonlinear regime. Note that our treatment can be applied to more general oscillating Joscphson circuit by 
simply renormalizing the value of A (for instance in the Cavity Bifurcation Amplifier (CBA)^SI or even an array of 
Josephson junctions inside a cavity^). 

Equating the currents at the node and expressing them in terms of the node flux we get the following equation 
of motion. It is the equation for a nonlinear, damped driven Duffing oscillator. Note that we have expanded the sine 
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FIG. 2: Characteristic frequencies of a double-pump amplifier (DPA). The upper/lower axes represents the frequencies at the 
input/output port. Diagonal dashed arrows represent the parametric coupling between various frequency components. The 
shade of the gradient in the middle band is a schematic representation of the oscillator bandwidth 2F. The resonant frequency 
of the oscillator is u>o- The pump frequencies are denoted by u>i and u>2 respectively and their average uj g defines the band center. 
The reduced detuning of u>o from u g is Q — (u>o — ui g )/T. The generic frequency of the signal and that of the corresponding idler 
are denoted by uis and uji — 2ui g — u>s, respectively. The reduced detuning of the signal (idler) from ui g is A = (uis — w g )/r 
(—A for the idler). The reduced detunings Q and A are not indicated to scale and have been exaggerated for clarity. 



function in JTJ current to third order for obtaining this equation 
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The source terms on the right are the drive currents - Irf^ Irf 2 , contributed by two pumps and the 'noise current' 
In, whose presence is imposed by the fluctuation-dissipation theorem. The damping constant T = \/2RC and the 
natural frequency of oscillation of the circuit is uq = \j\JhjC. Here A denotes the nonlinearity coefficient, whose 
value is —1/6 in the case of Josephson junction (cf. Fig. |3J. Also, the pump frequencies u)\ and u>2 are separated 
from ujq by several line widths T. 

To solve the above equation, we will use the formalism of input-output theory (IOT/221 where the electrical voltage 
and current signals at the level of resistance are treated as fields and are decomposed into incoming and outgoing 
traveling waves (cf. Fig [4]) propagating along a transmission line. The equations relating the voltage and current 
at any given position of the line to the wave amplitude (^4) of the wave propagating along the line are given by (cf. 
Appendix [b| : 



V(t) = V in (t) + V ont (t); 

p^in /out /~Z~ J^ in / ou ^- 



i(t) = i in (t) - 1 

in /out ]_ 



out (t) 

^in/out 



(2) 
(3) 



Here the wave amplitude A has dimensions of (Watt) 1 / 2 . We will invoke this technique to solve our final system of 
equations. 

Returning to the main Eq. (JlJ, we rewrite it as 



if + 2Tip + w 2 v(i + \<p 2 ) — w, 



2 iRFx 2 iRF-z 

Uq—t— 



(4) 
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FIG. 3: Top panel: Josephson nonlinear oscillator pumped with two current sources Irfx and Irf 2 - The Josephson element, 
denoted by a cross symbol, is used for its nonlinear inductance (SLj — (Lj/6)(L 2 J I 2 )/(h/2e) 2 + 0(I 2 )) which performs coherent 
multiwave mixing. Bottom panel represents full schematic of the DPA. The current Is corresponds to the input signal while 
In represents the noise current due to the resistance R, modeling the internal resistances of current sources. The symbol $ is 
the node flux associated with the Josephson element. 



where we have introduced a dimensionless variable ip = . The symbol v m = q^yf^ represents the quantum noise 
field driving the oscillator in the propagating mode picture. We assume a solution of the form 



ip- 



/oo 
X[u>] exp(—iuit)dui + {Sexp(— iui g t) + Y exp(— iui\i) + Zexp(— iu>2t) + c.c.} (5) 
-oo 



where 3 is a phasor corresponding to internal coordinate of amplifier oscillating at frequency u> g (= ^+" 2 ) . This 
describes the possible self-oscillation amplitude of the system. The symbols Y, Z are phasors corresponding to the 
pumps at frequencies uj\ and L02 respectively. Note that here we treat the generic signal (at frequency u)s) as a Dirac- 
delta component Xs[u;,s] = xs5{uj — los) + x* s S(uj + u>s) in the noise amplitude X[uj]. For further use, we introduce 
the "noise back-action factor" II as, 

/OO />OG 
/ duj a du b X{uj a )X{uj b }5(uj a + u b - 2u g ) (6) 
-00 J — 00 

Note that II is dimensionless like tp and unlike X[uj] which has dimensions of time. 
Substituting Eq. Q in Q, using the definition 



f(t) = J cLuf[uj] exp(-iujt) 

and performing harmonic balance (for frequencies wg, loj, u>x, L02) leads to the following system of coupled equations 
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FIG. 4: Traveling wave representation of circuit of Fig. [3] The resistance R is replaced with a transmission line of characteristic 
impedance Z c — R and the voltages and currents are treated as particular linear combinations of incoming and outgoing wave 
amplitudes. 
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Also note that X[wj] is the component of input noise at the idler frequency ujj 



with X = 
(uji + UJ 2 

As evident from the Eq. [6] the back-action factor II involves the mixing of signal and idler waves, realized through 
the nonlinearity of the system. It captures the effect of the reaction due to signal idler coupling on the pumps and 
eventually leads to pump depletion. Hence, it serves to make the coupling between pumps and signals symmetric in 
the sense that as the signal grows due to transfer of energy from pump into signal channel, it leads to a reduction 
of the pump amplitudes (as can be seen from Eqns. |9]and [l0| ) and this depleted pump then acts on the signal. The 
delta function in frequency in the expression for II is included to enforce harmonic balance (cf. Eqns. ^ and (10)). 

In deriving the above set of equations we have ignored the 2wo components and higher harmonics under the rotating 
wave approximation (RWA). Also note that we have neglected additional terms of the type |X| 2 + |y| 2 + |Z| 2 in each 
of the above equations as their effect on dynamics is a simple 'renormalization' of the plasma frequency (it will be seen 
later in section IIIB that u>o is actually a function of self-oscillation amplitude S and noise back-action factor II.) We 
have also dropped the terms of the form like IIX*[u;s], denoting explicit coupling between different back-action terms 
exclusive of the pump, as we will restrict our analysis to the perturbative limit where such terms are much smaller in 
magnitude. Under the above scheme of approximations, the system of nonlinear coupled Eqns. (6pTl form the basic 
framework of the problem to be solved. 



Pump amplitudes: Signature of back-action 

We solve Eqns. ([9| and (10 1 simultaneously, in the high-Q limit (r -C ujq), to obtain: 



yz = hh 



■1 + 7 n 



(12) 
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where 7 = , , , , U = -r^- /-, , 03 — ^tt are dimensionless constants of the problem. As is evident from 

the above equation, the correction to the pump amplitudes due to back-action (II and S 2 ) depends mainly on the 
nonlinearity parameter A. 



Eq. (12 1 gives the magnitude of the product YZ instead of individual pump amplitudes. It is convenient to 
formulate the drive strength in this manner as it is the product of two pump amplitudes which enters the equations 
of motion of signal and idler (cf. Eqns. Q, ([8])). Also, it should be noted that as we work in the perturbative limit 
or weak nonlinearity, we have included terms only up to the linear order in A in deriving the above equation. 

III. STIFF PUMP APPROXIMATION 
A. Signal and Idler Response 



Under the stiff pump approximation we ignore the effect of the back-action terms 7 (n + |^ 2 ) in Eq. (12). Thus, 

the pump amplitudes are simply given by YZ = — /1/2. This approximation will be further validated when we 
explicitly calculate the self-oscillation signal in subsequent sections. Using this value of drive amplitudes in Eqns. ^ 
and ^ and requesting to ~ ujq (under RWA), the signal and idler equations yield 

(n - A - i)X[u s ] - FX[-uj] = (J^j i/> s ] (13) 

(fl + A - i)X[ujj] - FX[-u s ] = (J^J (14) 

A = " s r " a represents the reduced detuning of the signal from the ghost frequency uj 9 . In this scheme of measuring 
the frequencies from the band-center uj g , the idler corresponds to a detuning —A from the ghost frequency. In the 
following expressions, we switch to the convention where = X[|cj|] and X[— = X*[\oj\] in order to avoid 

confusion between negative frequencies and negative reduced signal detuning. Also introduced is the reduced pump 
detuning il — "° r " g , and the effective pump power 

2A Irf\Ir.F2 

(for u>i — 3/4^0 ,^2 = 5/4cl>o)- On solving the two equations simultaneously, we get 

2 \ (O + A + i)v in (A) + Fv* in (-A) 



(15) 



X(A) 



X(-A) = 



uj q J (n-A-i)(n + A + i)-F 2 
2 

2 \ (fl - A + i)v in (-A) + Fv* iD (A) 



[A(A) v in (A) + B(A) v* in (-A)] (16) 



u J (n + A-i)(n-A + i)-F 2 
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w J 



) [A(-A) v in (-A) + B(-A) v* in (A)] 



(17) 



A plot of the |A(A)| 2 (cis-gain) and |i?(A)| 2 ( iraras-gain) , as a function of detuning from the ghost frequency, is 
shown in Figs. |5ja) and |5jh) for different pump powers F. This shows that the maximum gain is realized at zero 
detuning (at ghost frequency) and not at the natural frequency of the system w in presence of the drive. Internal 
gain profile at the signal frequency exhibits a frequency-pulling effect which comes into play as the drive is increased 
from zero. This leads to a progressive shift of the spectrum towards to g with increase in drive power and leads to 
maximum parametric gain at the bifurcation threshold. Note that the irons-gain does not 'see' the natural frequency 
ojq and maximal amplification occurs only at the ghost frequency for all drive powers. 

The coefficients -A(A) and -B(A) crucially decide the gain and noise performance of the amplifier. This follows from 
the relation between the input and the output fields of the amplifier 



v out [uj} = -iu(p[u] - v" 
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FIG. 5: (a) Signal-to-signal response of the internal coordinate to the incoming signal (internal cis-gain of the DPA) and (b) 
signal-to-idler response of the internal coordinate to the incoming signal (internal trans-gain of the DPA) plotted for reduced 
pump detuning Q = 2 as a function of reduced signal detuning A. The values A = 2 and A = correspond to u>s = and 
u>s =w g , respectively. The two spectra diverge in a similar fashion as the system is driven to the bifurcation threshold power 
F B . 



which leads to: 
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( a in [iv s } \ 
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V« in h^/] ) 



(18) 



where r — —2iA — 1 and s = —2iB. Thus, from Eqns. (16 1 and (17), we obtain 




(19) 



Here we have introduced a ln / out are the respective field operators for the input and the output ports of the line (cf. 
Appendix |b|) . It can be verified that the resulting scattering matrix is symplectic in nature, i.e. T SJS = J2H where 

/ 1 
-10 
1 
\ 0-10 

and has a unit determinant. The property of symplecticity ensures the absence of any extraneous degrees of free- 
dom, an d hen ce no missing information (information preservation) - a condition necessary for quantum-limited 
detection 29 30 . Moreover, Irl 2 — \s\ 2 = 1. Thus we can write 



a out tas] 



G a in ta s ] + y/G^l a in [-U!} 



(20) 



when ujs and luj are very close and phase factors can be ignored. Here G — \r\ 2 is the power gain of the amplifier. A 



relation of the form shown in Eq. (20 1 is typical of an amplifier operating at the quantum limit. The internal mode 



fluctuations of the amplifier are the vacuum fluctuations at the idler port, which are the source of noise added by the 
amplifier. The efficiency of their conversion to signal frequency at the output port is indicated by the inter-conversion 
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FIG. 6: Upper panels show gain profiles for signal (|r| 2 ) and inter-conversion (|s| 2 ) gains while lower panels show cis (\A\ 2 )- 
and trans (|B| 2 )- gains as a function of reduced drive strength F and signal detuning A. Here the value Q = 5, which is 
the upper limit of the frequency range in Fig. [5] was used to display the salient features clearly. The plots for \r\ 2 and |s| 2 
are symmetric about the ghost frequency corresponding to A = 0. The internal trans-gain of the amplifier also exhibits this 
symmetry. However, the cis-gain is asymmetric since it knows about the relative positions of ujq and u) g . 



gain (or |s| 2 ). Fig. [6] shows the profile of the signal and inter-conversion gains on a two dimensional color plot, as 
a function of drive power and detuning. For the sake of comparison, plots for cis- and trans-gains are also shown. 
The maximum signal and inter-conversion gain profiles are symmetric and indicate maximum gain for both positive 
(signal) and negative (idler) values of detuning at the ghost frequency and the bifurcation threshold. 
Further, we can calculate the output noise spectrum from Eq. (20 1 as 



-rin/out 

Hlo s {G- l)SZYoi] with S'Z = icoth(^) (cf. Appendix 
noise temperature of the amplifier, T/v 



where 5 o / ° u '' represent the input/output photon number spectral densities. This gives the added noise as 

B|. When referring back to the input, we obtain the 



fe B G' 



as 



T N = 



2k h 



(21) 



for G>1 and T — > 0. Thus, like a single pump paramp, the DPA adds half a photon at the signal frequency and 
approaches arbitrarily close to quantum limited behavior when operated at ksT <C ftw. 
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B. Steady State Calculation 



Sclf-oscillation is defined as the response of the system in the absence of any input. To compute this for the DPA, we 
use Eq. ( 12 ), under the stiff pump approximation (ignoring the noise back-action factor II). We need to reinstate the 
E! 2 term in F for this calculation as the coupling between 3 and F is responsible for development of a self-oscillation 
amplitude in the system. Thus a strictly stiff pump, like that used for derivation of noise spectra in the preceding 
section, cannot lead to spontaneous emission of finite amplitude. The inclusion of back-action to some zeroth order, 
i.e. due to self-oscillation itself and not due to noise, is crucial for the prediction of this effect. 

From Eq. (Ill, we have 



(n-i)z + z*F(-i + ^z 2 ) =o. 



(22) 



Writing the equation for complex conjugate and eliminating s* from both the equations in favor of s, we get a 
cubic equation in |S| 2 . 



1 



(fl d + 1 - F Z )\E\' + ft 7 F|S| 4 + -YF"\S\° = 



(23) 



On factoring out the zero solution |S| 2 and ignoring the imaginary solution, we get the solution for self-oscillation 
amplitude as: 




(24) 



A plot of Eq. (24) (cf. Figs. [7^i and [7)3) shows that beyond the threshold power Fg, the system develops a finite 
amplitude of sclf-oscillation. This phenomena is reminiscent of lasing although here the frequency of self-oscillation 
is imposed by a combination of the frequencies of the two pumps. This also justifies our simplifying assumption of 
ignoring the back-action due to S while calculating signal and idler gains, as we see that before bifurcation threshold 
this approximation is exact. 

The phase diagram of the DPA (cf. Fig. [8^) depicting the locus of bifurcation follows from Eq. (24 1. Since 
F B = Vft 2 + 1, we see that the phase diagram is symmetric about 17 = 0. Also, the minimum pump power required 
to make the system self-oscillating (F c ) is unity, which is realized when system is pumped at resonance (wi+cja = 2o>o)- 
This scenario may be contrasted with the phase diagram of the Josephson bifurcation amplifier (JBA) which employs 
a single pump source (cf. Fig. ^p). The phase diagram of the JBA shows an inherent skewness, which is absent in 
the DPA. Moreover, the optimal bias point for self-oscillations does not coincide with the optimal point of operation 
ensuring maximum parametric amplification. For the DPA, by contrast, the two points coincide at the global minimum 
of the F — plane. Appendix [X] contains a detailed discussion of the JBA phase diagram. Calculation of explicit 
transition probabilities between the two stable states resulting in amplification near bifurcation for single pump 
systems can be found in Ref. 31. 

It is useful to plot the linear response of the amplifier to wide band noise and classical self-oscillation together, as 
in Fig. [9] The integrated noise power gain is calculated as 



P out 



JV 



1 

2T 



\r[v]\*d(w). 



(25) 
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Fig. [9]shows an important result: the threshold power for lasing coincides with the power corresponding to maximum 
of the integrated noise spectrum. This strongly suggests that the origin of self-oscillation of the system is a result 
of coherent multi-wave mixing of correlated noise components, which grow around the effective pump frequency uj g 
with an increase in the drive. These noise correlations, resulting from parametric interactions, grow stronger as we 
drive the system near bifurcation threshold and ultimately lead the system to self-oscillate. Thus, beyond bifurcation 
threshold, there is an additional channel coinciding with the effective pump frequency oj g in which the system leaks 
energy. 



C. Squeezing 



One of the most a ttractive applications of parametric amplifiers is the generation of squeezed states of electro- 
magnetic signalp2122l A squeezed field is one with phase-sensitive quantum fluctuations, which are amplified in one 




FIG. 7: Self-oscillation amplitude as a function of drive F and reduced pump detuning Q. Panel (a) represents cuts of the 
surface shown in panel (b). Note that self-oscillation develops only beyond a bifurcation threshold Fb which has a parabola- like 
dependence on the reduced pump detuning. The lowest bifurcation value Fb = 1 is found at zero pump detuning (i.e. uj g = Uq) 
and defines the minimum/critical bifurcation power F c . 
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FIG. 8: (a) Phase diagram of the DPA in absence of noise at the input port. Self-oscillation occurs above the solid line. The 
vertical dashed line represents the condition for optimal parametric pumping. In this diagram, the bifurcation value Fb is 
a single- valued function of Q. (b) Corresponding phase diagram of the Josephson Bifurcation Amplifier (JBA), showing the 
bifurcation curves for upper (blue) and lower (pink) bifurcation points. Here, there exists a critical detuning S7 C = V3 below 
which the system does not exhibit bifurcation. Below f2 c , the bifurcation power is a double-valued function of f2. The curved 
dashed line corresponds to the location of maximum oscillation amplitude. See Fig. 18 (Appendix |A| for the behavior of 
oscillation signal as a function of drive strength in the JBA. 



quadrature and deamplificd in the other below the quantum vacuum floor. They exhibit strong nonclassical effects 
which find applications in optical communication systems, interferometry-based gravitational wave detection^! and 
quantum cryptog raphy 3 -^. S queezing using Josephson parametric devices has already been realized in four-wave mix- 
ing configurations^*- 37 * 38 * 39 ! using single pump. Here we will show the identification of squeezing produced by the 
two-pump scheme. 

The basic idea is to exploit the correlation between the generated signal and idler photons. The output field of the 
oscillator is first mixed with a signal at ghost frequencj®]. Then, the homodyne spectrum at the output of the mixer 
can exhibit sub-quantum noise level for a suitably chosen phase difference (6) between the original drive and beat 
signal. Using the standard quantum optics procedure, the mixing operation in frequency domain can be written as: 

a mix (A) = exp(-i(9)a out (A) + exp(i(9)a out (-A). (26) 

It can be confirmed that the above relation is identical to that obtained for t he o utput field operator of a two-mode 



squeezed vacuum state, under the identification G(— \r\ ) — cosh t in Eq. (20 PM Using the expressions for a' 
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FIG. 9: Integrated noise power gain and self-oscillation amplitude, for the case of stiff pump with = 1.0, as a function of 
reduced pump power. The pump frequencies were uj — 3/4cjo,^2 = 5/4wq, giving 7 ~ 0.1. 



coupled with the commutation relations (cf. Appendix |B|) , 

({a in (A),a^(A')}) = S aa S(A-A') 
where S aa is the spectral density of quantum noise level, we get 



- m,x - S* aa (|A out (A)| 2 + |B out (A)| 2 + |A out (-A)| 2 + |£ out (-A)| 2 

2i?e{exp(2i6»)(A out (A)B out (-A) + A out (-A)B out (A))}) (27) 



a a 



Here A out = -2iA ~ 1 and B out = -2iB with A and B defined by Eqns. ([16J and ^7j. Thus, S™* is reduced below 
Saa for a particular value of 9 and shows a maximum for the orthogonal quadrature. The squeezing spectra for the 
amplified and deamplified quadratures are shown in Fig. [10] 

It can be seen that maximum squeezing is attained at zero detuning (lus = w g ) and the bifurcation threshold 
corresponding to a given It may be noted that perfect squeezing at bifurcation is obtained as a result of our stiff 
pump analysis. The more sophisticated soft pump analysis would predict a squeezing fraction less than 100% even at 
the bifurcation. 



IV. SOFT PUMP: INCLUSION OF BACK-ACTION CORRECTIONS 



To incorporate the effect of back-action on the dynamics, we use the full expression given in Eq. (12 1 for the 
pump amplitudes to calculate the response functions of interest. We adopt a mean-field approach in this calculation. 
Various levels of coupling in the system, manifest in Eqns. ^ to (111, necessitate such a self-consistent analysis of 
the problem at hand. To make a reasonable start, we pick up one of the strings of this self-consistent loop - signal 
and idler amplitudes, and calculate them treating LI as a parameter of the problem to be evaluated later. 
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FIG. 10: Squeezing action of the DPA. ( Upper panel) Deamplified and amplified quadrature power gain for different drive 
strengths, as a function of reduced signal detuning, with Q — 1. Maximum squeezing corresponds to the bifurcation threshold 
(F — Fb in Fig. [9|, obtained for both peak noise power and onset of self-oscillation. (Lower panel) Squeezing fraction = 

1 — ^p^j increases from (no squeezing) at zero drive to its maximum value 1 (perfect squeezing) at ^-dependent bifurcation 

threshold. Beyond the threshold, squeezing decreases. 



A. Corrections to signal and idler amplitudes 



Using the same method as in the previous section - this time using the full expression of pump amplitudes from 
Eq.(12l - we solve for signal and idler amplitudes and obtain 



X[w s ] = 



2_\ (fl + A + i)v ia [ujs} + F (1 - 7 (n + j^S 2 )) v^I-lj!] 

J (fi — A — + A + i) - f 2 (i - 7 (n + ^s 2 )) (i - 7 (n* + |s* 2 )) 

2 



CJ 



A[w 8 ] v m [u> s ]+B[u;s] ^[-wj] 



(28a) 



X[cjj] = 



(n-A + iV n [u;j] + F (1 - 7 (n + |S 2 )) v in [-m s } 

(n + a - - a + i) - f 2 (i - 7 (n + is 2 )) (i - 7 (n* + is* 2 )) 



(28b) 



where A and B are the corresponding coefficients for a soft pump, analogous to those in Eqns. (16 1 and (17 1. It is 



convenient, both in terms of notation and analysis, to define a back-action corrected drive power, Fba, given by 



F, 



BA 



F 



1-7 n 



1 



— 2 



(29) 



Note that Fba, unlike F, is complex and depends self-consistently on X[ui] through an integral (II). 

As before, we calculate the signal and inter-conversion gains using u out [w] = — iwA"[w] — t> m [a>] in terms of this new 



drive strength. The resulting scattering matrix S, defined by a° 



Sa ln [u>] is again symplectic despite the complex 



nature of Fba- This shows that the basic tenet of information preservation requested during an amplification process 
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is fulfilled even in the presence of back-action. The scattering matrix of the amplifier is evaluated to be 



S 



( e ia r 

e~ ia r* e ~ i ( a + s ) 

e- i(a+s h e~ ia r 



s i(a+<5) , 






\ 



with 







(q+(5) s * q 





e ia * J 

/ 




-n 2 


+ A 2 + |F SA | 2 ] 






w 


- A - 


+ A + i) - 
~2iF BA 






s m 


- A- 


+ A + i) - 






a — arg 




1 




; 8 = wg[F BA \. 


(0 - 


A -i)(Q + A + 


i) - I^baI 2 . 



(30) 



As before (cf. Section III A), \r 



1. Note that the expressions in (30) reduce to those of Eq. (19 1 when 



Fba — F. The presence of the additional phase factor S, which is a function of pump detuning f2, can be understood 
by noting that II is a complex quantity in general. 



B. Corrections to Steady State Signal 



We recalculate the steady state response (i.e. zero input response) of the system, using the back-action corrected 
drive strength. Eq. (Ill gives 



(n-i)~ + x ~*F 



-i + 7 n 



i 



— 2 



0. 



On solving this equation, as in section ( III B I , we get the self-oscillation amplitude as: 



|S| 2 = 2 



-M + ^\F BA \ 2 -\ 



(31) 



(32) 



It is instructive to see that Eq. (32 1 reduces to Eq. (24 1, in absence of back-action i.e. Fba 1 — * F 



V. CALCULATION OF THE NOISE BACK-ACTION FACTOR II 



We will use the results derived in section IV A to calculate the noise back-action factor II. Using Eqns. (28a I and 
( 28b I in Eq. ^ , under mean field 



n 



OO pQO 

duj, 

2 



duj b X[u) a }X[uj b ]5(uj a + u) b - 2u> g ) 



du> a / d,LO b {A[uj a ) v m [uj a ) + B[ui a ] v m [uj a - 2uj g )}x 



{A[u b ] v in [uj b } + B[w b ] v in [uo b - 2uOg]}5{u a + uj b - 2uj g )) 



Exploiting the fact that (U m [cL> a ] • V 111 ^]) = S vv [uj a ]8(ui a +u b ) (cf. Appendix 
two terms out of the above four survive giving 



B| and recalling that v ln — ^ 27r » on ^ 



/ 4tt \ 2 f°° 

n=[^-J J duj a {A[uo a ]B[u b ] + B[uj a ]A[uj b ]}S vv [uj a ]. 



(33) 



Using 



Z c haj a 



coth 



2k B T 



T »w^ Z c k B T 
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in Eq. (33), we obtain 



/oo 
d(A)V(Sl,A,F) 
-oo 



with 



V(Sl,A,F) 



2(si + i)F BA 



e 



eff 



{(SI - A + i)(Sl + A - i) - F BA }{(Sl - A - + A + i) - F| A } 
k B T 
Ej ' 



(34a) 

(34b) 
(34c) 



Here we have used E 



J 



r, r 



i 

IRC' 



To evaluate II in Eq. (34a), we use the method of residues. We find that the denominator of V has four poles in 



the complex plane of reduced frequency A at (— i ± y^Sl 2 — Fg A ) and corresponding complex conjugates - giving one 
pole in each quadrant (Fig. 111. To respect causality we need to take into account, at any given drive power, only 
the contribution of the poles in the upper half complex plane. We evaluate the residue at each of the relevant poles 



and sum them up to get the value of complex- valued integral in Eq. ( 34a I and obtain 
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Re A 
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FIG. 11: Poles of the noise back-action factor II in the complex plane of reduced signal frequency A. They move as the drive 
strength is changed. The poles in panel (a) before the bifurcation threshold move towards the imaginary axis. In panel (b), 
they hit the imaginary axis as Fba takes the value Q.. As Fba is further increased, in panel (c), they move along imaginary axis 
in two opposite directions. Finally, in panel (d), one pair of poles hits the real axis as Fba attains the bifurcation threshold. 



//7T 



n = 



?/7T 



[SI + i)F BA 

1-F% A + Sl* 
(SI + i)F BA 



("1 + ^ 



Sl*)jF? 



BA 



SI 2 



f ba < VsWTT 
f ba > VsWTi 



(35) 



The expressions for Id are continuous at the bifurcation threshold but have different slopes and qualitative de- 
pendence on Fba- This indicates that the response functions calculated using II, though continuous may not be 
symmetric across threshold power as in the case of stiff pump. Also introduced is a new parameter 



rj = 7 • 6 e ff 



(36) 
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which occurs as a pre-factor of II in various gain coefficients and incorporates the effect of both non-linearity(7 oc A) 
and effective temperature of the input (0 c ff oc T) . This serves to restrict the parameter space of the problem and can 
be identified as a kind of back-action index. Henceforth, we will define the system parameters in terms of r\. 

As evident from Eq. (35), II depends on Fba, which itself depends on II. Such a relation is a natural consequence 



of a mean-field approach. Also, it should be noted that as the poles depend on this corrected value of drive power. We 
need to fully appreciate the three-way coupling reflected by Eqns. (29), (32) and (35) before embarking on a solution. 



Both the complexity and strength of the following analysis lies in the self-consistent treatment encompassing all the 
three players of the game. In these three equations, we identify Fba as a convenient "slave" parameter (Fig. 12) and 
solve the system of equations parametrically. 




Eq. (35) 



FIG. 12: Schematic diagram displaying the structure of the set of Eqns. ( 29 \ , ( 32 l and (35 I. The drive strength F, self-oscillation 



amplitude 3 and noise back-action factor II, all involve back-action corrected pump strength Fba and are interdependent 



It may be noted that the self-oscillation term is especially important once the bifurcation threshold is reached and 
needs to be incorporated to get the correct value of bac k-action on drive strength. However in the expressions for 



signal and inter-conversion gains - \r\ and \s\ (Eq. 30 1- only terms like I-FbaI appear. Due to this, only the real 
and absolute values of n and S 2 are important. Exploiting this fact, we have used the expression obtained for the 



self-oscillation signal as obtained in Eq. (32 1, to account for the back-action of self oscillation signal on drive strength 
(S 2 term in Eq. rfggt). 



In Fig. 1 1 3 1 we show parametric plots of both the integrated noise power gain and self-oscillation amplitude of 
the system, obtained after incorporating the effect of back-action. The response functions clearly show that the 




FIG. 13: Maximum signal gain (dotted curve), integrated noise power gain (solid curve), self-oscillation component (dotted- 
dashed curve) as a function of reduced pump power for f2 = 0.5 and n — 3 x 10~ 3 (values corresponding to |A| = 1/6 and 
@eff — 3%). The dashed curves show the respective stiff pump responses. Back-action manifests itself through a pronounced 
reentrant behavior and a shift of the bifurcation threshold to higher powers. The inset shows the noise back-action factor II 
as a function of F evaluated for the same system parameters. The slope represents the impedance seen looking into the pump 
port. This impedance acquires power dependence as bifurcation threshold is approached. 
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FIG. 14: Integrated noise power gain in dB (solid curve) and self-oscillation (dotted-dashed curve) responses, at Q = 0.5, for 
different values of temperature (a) r\ = 10 -3 (b) 77 = 3 x 10 -3 and (c) r\ = 10~ 2 . The noise back-action factor scales linearly 
with temperature of the noise at the input and hence its effect becomes more pronounced at higher temperatures. Curves above 
reflect this through a higher fold-over and increase in threshold power with temperature. The dashed curves representing the 
behavior of respective stiff pump response functions are included for comparison. 




FIG. 15: This plot shows the dependence of system response on pump detuning f2 in the presence of back-action for r\ = 10 . 
Back-action manifests itself in the form of hysteresis whose range increases as we move away from a renormalized resonant 
frequency of the oscillator, whose value corresponds here to Qq = —0.03. This value, at which hysteresis disappears, is offset 
from the stiff pump lowest bifurcation value Q — 0.0. As in Fig. |14| the dashed curves representing the respective stiff pump 
behavior are included for comparison. 



system exhibits 'hysteresis' as indicated by the reentrant gain curves: gain is a multiple valued function of the control 
parameter (drive strength here). Another consequence of the incorporation of back-action is the reduction of the 
effective drive power which shifts the threshold power corresponding to the maximum noise gain, as well as the onset 
of self-oscillation, to higher values. Thus, the behavior of system in the presence of back-action is markedly different 
from that obtained in the stiff pump approximation (Fig. KjJ. 

We now study the effect of various system parameters on the back-action induced corrections. An increase in 
effective temperature is seen to cause a greater back-action correction. This is depicted in the noise gain and self- 
oscillation gain curves plotted in Fig. [14] for different values of parameter 77, which varies linearly with effective 
temperature of the input port ksT/Ej (cf. Eq. ( 36 1 ) . 

Another parameter of interest in the problem is the reduced pump detuning Q. ft is seen that back-action pushes 
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FIG. 16: Variation of hysteresis due to back-action in integrated noise power gain as a function of pump frequency. The 
direction of back-bending of the response curves reverses as we move from negative to positive detuning. 



the bifurcation threshold powers to higher values as f2 is increased or as we move away from the resonant frequency 
wq- Besides an increase in threshold power, the extent of fold-over of gain curves as a function of drive strength 



also varies as we vary (Fig. 16 1. This is depicted in Fig. 15 An important ramification of this effect is that 
£1 = 0(iUg = ujq) does not correspond to resonance any longer, as in the case of stiff pump. This is a consequence 
of the fact that the resonant frequency of the system lu (defined as \/\JLjCj ) itself changes in the presence of a 
drive, due to back-action. This indicates that for some |f2o| — 0.03, the system simulates stiff pump behavior and the 
corresponding luq = uj g + CIqT corresponds to a corrected resonance frequency. Thus, by tuning f2 close to the new 
resonance, hysteresis can be eliminated. However, to realize maximum parametric gain at optimal powers and relieve 
the system of hysteresis we need to tune both the parameters F and f2 to ensure the operation of the amplifier near 
the optimal point in the phase space. 

This can be further affirmed by looking at the phase diagram of the system. The locus of points showing the new 



bifurcation thresholds as a function of fl determine the new phase boundary (Fig. 17). We note a shift of the phase 
boundary on incorporation of back-action which causes the minimum of the phase-transition curve to move away 
from Q = (Fig. [8| towards negative frequency axis. This occurrence can be attributed to the decrease in resonance 
frequency loq with an increase in amplitude of oscillation due to higher drive strength (cf . Eq. ( A3 ) for the typical 
form of dependence). This causes the effective detuning — (luq — u> g )/T assume a non-zero negative value in the 
presence of noise. The above shift also explains the switch between the two branches of the bifurcation loci shown in 
Fig. 17 for a given r\. The system traces the upper (thin) curve showing a decrease in oscillation with drive, till the 
point it traverses through the minimum where it become unstable and switches to the other branch of high amplitude 
oscillation (thin portion of the other curve) where the resonance condition can be better met. 

Another important result highlighted by the phase diagram is the necessity of biasing the system near resonance. 
A comparison with the curves plotted in Fig. [14] shows that the system is much more robust to back-action if biased 
near the relevant frequency (new loq corresponding to the shift in il). For instance, c ff = 0.1 corresponded to a large 



shift in threshold power for bifurcation ~ 20% (cf. Fig. [Tip ) while the corresponding curve in Fig. 17 shows a shift 
of 5-10 % in bifurcation threshold, if biased near new resonance. 

It may be noted that for higher values of f2(> 1), the frustration in the system increases further and the resultant 
corrections arising due to back-action approach order unity leading to anomalous behavior. Hence the underlying 
assumptions of perturbative treatment of back-action demands that we restrict ourselves to cases where the system 
is pumped sufficiently close to plasma frequency of the junction (\luq — Lu g \ <C 1). 



VI. CONCLUSIONS AND PERSPECTIVES 



We performed a first-principle analysis of a Josephson parametric amplifier involving two RF energy sources and 
calculated the reciprocal effect of pump amplitude on signal and idler gains in a self-consistent manner. The analysis 
was based on a mean field approach of the intrinsic couplings between various components of the frequency spectrum, 
especially near bifurcation threshold. In particular, the effect of back-action on integrated noise power gain and 
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FIG. 17: Phase diagram for the DPA with back-action correction, for different values of effective temperature Q c g. At each 
temperature, the curves corresponding to the two branches of noise back-action factor II are plotted parametrically (cf. Eq 35 1 . 
The region between the two curves corresponds to the hysteresis of the system. The intersection of the two curves marks the 
locus of second-order phase transition points (dashed line). The bifurcation threshold is seen to shift upwards to higher pump 
strengths and towards higher values of Q with temperature. 



sclf-oscillation amplitude response of the system were investigated. 

Starting with the equations derived in sect ion |IT| response functions for the DPA were calculated in section |III| 
for a stiff pump. Following analysis in section |III A| showed that the DPA can indeed achieve quantum limited noise 
temperature. The DPA also exhibited squeezing with fluctuations in one of the two quadratures squeezed below 
the quantum limit (cf . section III C I . Maximum squeezing was found at the bifurcation threshold for a given pump 
detuning from resonance, as expected. Then in section [PV] corrections to signal, idler and self-oscillation signal gains 
were calculated, in the presence of full back-action, yielding our main results as set of coupled equations - (29 1, (32 1 
and (35 1- involving three inter-dependent response functions of the system S, II and signal, idler amplitdes. 

The significance of the inclusion of back-action is manifested even at the stiff pump level, as seen in section |III B| 
for the calculation of sclf-oscillation amplitude. It was shown that the inclusion of back-action to zeroth order (i.e. 
back-action due to self-oscillation signal only) is imperative to get steady state response of the system. These results 
are relevant for the amplifier designs based on three-wave mixing with SQUIDS as the inherent flux-based coupling 
mechanism impels a self-consistent treatment of pumps while calculating signal gains. We believe that the recent 
experimental results showing a spontaneous emission, reported for a flux-driven Josephson parametric amplifiei!^, can 
be well explained by doing a similar treatment. Another important feature that emerged out of this restricted analysis, 
with a stiff pump, was the coincidence of maximum-gain and bifurcation threshold. This harmony is contingent to 
symmetric bifurcation boundary obtained for the DPA. 

A self-consistent treatment of back-action led us to observe hysteresis in the system, as a function of drive strength. 
The magnitude of hysteresis increases as the detuning of the ghost frequency lo q from resonance and effective tempera- 
ture at the input port of the oscillator are increased. However, the feature of coincidence of maximum parametric gain 
and self-oscillation thresholds is preserved (as in the case of a stiff pump). An important conclusion of our work is the 
shift of the resonant frequency of the system itself. This can be clearly deduced from the shift of the phase boundary 
of the system, which shifts to higher drive strengths and higher resonant frequencies with increase in back-action. 
This indicates that adjusting the effective pump frequency uo g near this 'new' resonance is crucial for making the 
system less susceptible to back-action induced hysteresis and realize bifurcation at optimal pump strengths. Future 
experiments with devices employing the bifurcation of Josephson junction should take this effect into account. Also, 
the fact that DPA shows quantum limited behavior in the limit of large gains can be exploited to build a traveling 
wave amplifier using a cascaded chain of single stage amplifiers, with the effective noise temperature of each stage 
reduced to at least the cooling chamber of the next stage. 

To summarize, we have developed a minimal model inclusive of all the various components required to understand 
the dynamics of microwave parametric amplifiers based on purely dispersive elements like Josephson tunnel junctions, 
both away from and near the bifurcation threshold. The power and usefulness of this analysis lies in its generality 
which makes the techniques developed for the analysis of the DPA readily applicable to a host of systems operating 
either as three-wave or four-wave mixing devices. 
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APPENDIX A: PHASE DIAGRAM FOR THE JOSEPHSON BIFURCATION AMPLIFIER (JBA) 



The equation of a damped RF-driven JTJ with a single pump, as in the model for the Josephson Bifurcation 
Amplifier (JBA), can be written as: 



d 2 d , 



Irf 
~C~ 



cos (tot) 



where A = 1/6, T = and — jjj— . As before, we assume a harmonic solution to the above equation 



(Al) 



(A2) 



which on performing harmonic balance and making the rotating wave approximation with 5(— uiq — to), T <C luq, leads 
to 



-ofi + 25lo + 2iLu T + Uq 1 - % |S| 



I , — . 1 2 



with x = x an< ^ f = ^oi^ 1 - ^ e can wr ite the above equation as 

S — aujQ |S| 2 ) + iT S 



J_ 
2uj 



(A3) 



where a = | is the anharmonicity parameter (variation of resonant frequency with energy) . Its value is 1 /SEj for 
the quadratic+quartic potential resulting from the expansion of the Josephson cosine potential 



uj rp = luq ( 1 — a IS 



O 



I , — .i 2 



(A4) 



Note that in the "Transmon" limitP^, this equation can be written as 

^12 = wqi (1 - atkj 01 ) . 
Using lu 01 = y/8EjEc /ti, we recover the useful result 

UJl2 — UJqi — Ec/h. 

We can further reduce the equation for 5 by introducing 

n=- f= QCJ °I S I 2 f= /2 a = 



a 



Using the above parameters, Eq. ( A3 1 assumes the form 



(ft -E) 2 + l 



E = F 



For comparison, the corresponding state equation for the DPA is shown below (cf. Eq. 23 1 

[4(fi 2 + 1 - F' 2 ) + 4flFE + F 2 E 2 ] E = 



(A5) 



(A6) 



(A7) 



(A8) 
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FIG. 18: Comparison between the steady state oscillation amplitude of the DPA (left panel) and the JBA (right panel), at the 
ghost and pump frequency respectively, in absence of noise. In both the case, these frequencies correspond to the average of 
the signal and the idler frequencies. Unlike the DPA, JBA displays a steady state oscillatory signal as soon as the pump is 
turned on. Note that the dimensionless pump power F is related to the physical pump strength by different numerical factors 
in the two systems (see Eqns. 15 and A6l. 



where E = j\^\ 2 - A plot of the steady state behavior of the JBA is shown in Fig. 18 The analogous plot showing the 



steady state response of the DPA (as obtained in section III F3 1 is also displayed for comparison in the same figure. 
The plot for the DPA shows that before the bifurcation threshold, there is no response unlike the JBA. Also, note 
that there is no hysteresis in the DPA unlike the JBA, for which there is a lower and an upper bifurcation threshold. 
Thus the nature of transition between dynamical states at bifurcation is markedly different for the two systems. For 
fixed F the reduced detuning f2 can be obtained as a function of the reduced energy E 



<> = e± 



(A9) 



Therefore (for real values of detuning f2), we dem and E, 



F. We can find the location of the bifurcation for a 



given F by requesting 4§ = 0. Differentiating Eq. JA7b, we get a condition for the extremum points 



dn [2E (n - e)] + dE -2E(n-E) + i + {n-E) 2 =o 



(A10) 



This shows that there exists a critical detuning and drive such that E as a function of B and f2 has a triple real root. 



This is the location of bifurcation in E — ft plane. This can be found by locating the extremum of Eq. ( A7l by using 



the above condition as a function of E and requesting that they coincide. This gives us an equation for E 

3E 2 - 4En + 1 + n 2 = 

with two roots, whose values are given by: 



E = 



2f2±^4Q 2 -3(l + ft 2 ) 
3 

2f] ± V^ 2 - 3 



They are degenerate for the critical values 



3%/3' 



(AH) 



(A12) 



In semiclassical terms, this critical point corresponds to the drive being such that the average energy leads to a 
detuning comparable to the line width. 

We reconsider Eq. (Al I, but now include a weak 'signal term' in addition to the pump 



d 2 d 

X + 2T — X + uj 2 X (1 - X X 2 ) = / cos (uj d t) + e cos (u) s t) 



dt 2 



(A13) 
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FIG. 19: Phase diagram of the JBA in E-fi plane, showing the bifurcation curves (dotted-blue), the curve corresponding to the 
most efficient parametric amplification (solid-green with dashed portion indicating the bifurcation regime), and the asymptote 
to upper bifurcation level(dotted-dashed, red). Note that the Q = 2F line does not intersect the bifurcation curve at the lowest 
value of detuning 0, = \/3 marking the commencement of bifurcation regime. 



with e <C /. We find a solution of the form 



X(t)=S (t) + SX (t) 



where, SX (t) = xcos [(lo^ + uj m ) t + <p x ] + y sin [(coj — aj m ) t + <f> y ] , aj m = u> s — <C oJd, SX (t) < S (i) . Under the 
above approximations, we can write the equation for SX as 



d 2 



5X + 2T — SX + uj 2 SX -3 x ^SX[E(t)} 2 = ecos(oj s t) 



d 2 d 

—SX + 2T—SX + ujISX - Zx^lSX 

(Jib (XT/ 



|S| 2 E 2 e- 2luJdt + c.c. 



e cos (u) 8 t) 



(A14) 



where in the second step, we have used Eq. (A2). Introducing 



= cjo ( 1 — 2a |H| 1 =^ lu 



we obtain 



dt 2 



SX 



2T— SX + lo 2 . s [1 + 2 |e| sin (2uj d + ip)] SX = ecos (uj s t) . 



(A15) 



This is just the equation for a parametric amplifier. Therefore, the line f2 = 2E denotes the location of the optimal 
drive frequency for maximum parametric amplification. Note that it is not the same line as the zero-phase shift drive 
frequency (fi = E, the upper asymptote of the bistability line) , nor the line passing through the location of the critical 
point (ft = \E), nor the lower asymptote of the bistability line (ft = 3E). In the DPA, the optimal drive frequency for 
maximum parametric amplification remains O = below threshold. The above bifurcation behavior is quite generic 
in nonlinear oscillators. The theory of these systems is fairly well developed and a detailed exposition of these ideas 
can be found in Ref. 43. 



APPENDIX B: INPUT-OUTPUT THEORY 



Input-output theory (IOT) is a particular model of scattering theory (S-matrix theory) which applies to a system 
coupled to a heat bath. It is well documented in the literature but we include a brief description here for the help of 
the readers and consistency of notation. For the analysis via IOT, the resistance (R) of a circuit is replaced with a 
transmission line of characteristic impedance Z c (= R) and the voltage and the current along the line are expressed 



in terms of superposition of incoming and outgoing waves (Fig. 20 ) . The waves represent either a signal launched on 
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FIG. 20: Top panel shows a damped LC oscillator driven by an RF current source. Bottom panel shows the equivalent circuit 
in which the current source and its internal resistance have been replaced by a semi-infinite transmission line (input output 
analog). The two cases are equivalent from the point of view of the LC oscillator if we make the identification Zc ~ R and 
A in = s/RI h /2. 



the line to drive the oscillator (pumps, signal) or the thermal/quantum fluctuations in the line (e.g. Nyquist noise of 
the resistor). The power of this semiclassical technique, apart from its calculational advantage, lies in the provision 
of simple physical insights into the link between the noise sources and dissipation. The voltages (V) and currents (I) 
are expressed in terms of incoming and outgoing field amplitudes (A) are expressed as: 



V(z, t) = ^JZ c (A out (z, t) + A in (z, t)) 



(A out (z,t)~ A in (z,t)) 



(Bl) 



It is straightforward to obtain constitutive relations linking input, output and internal fields of the amplifier by 
imposing the appropriate boundary conditions at the termination of the line (z=0) 



V(t) = V in (t) + V out (t) ; I(t) = I in (t) - j 

where we have used the relation V in ^ out = y/Z~ c A in / out and I in / out = 
It is useful to define the quantities a[oj\ as 



= A[w] 



(B2) 



(B3) 



where A[ui] = — ^= J dtA(t) exp(iwi) . This leads us to a natural generalization to quantum regime as the normalized 
field amplitude a plays the role of bosonic field operator as defined for a harmonic oscillator. It obeys the following 
well-known field theoretical commutation relation 44 



[a in [oj},a in [Lu'}} = sgn[w] 5(u + uo') 



(B4) 



The fluctuations of field creation and annihilation operators is characterized by the corresponding noise spectrum in 
thermal equilibrium 



({a in M, a h V]})T = Saa[w] 8{u> + J) 

1 1 



S„ 



l , h\uj\ 



(B5) 
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The quantity SwS 1 ™ denotes the total energy per mode and reduces to ^f- in the limit of zero temperature (vacuum 
fluctuations) and the classical limit of k B T in the limit of high temperature. Eqns. (B4| and (B5l are valid over the 
entire frequency range, including the negative frequencies. We can return to the conventional description restricted 
to only positive frequencies by the identification 

a[— uj] — > a^[w]. 

The preceding equations lead us to define the ordered spectral density 



tlU) 



coth 



\2k B T 



5{uj + u)'). 



(B6) 



Thus, we can easily write the fluctuations of the voltage across the resistor 



<VMVV]) = 



z r hu 



coth 



huj 
2k B T 



1 



5(u} + u/) 



(B7) 



which follows from y in /° ut = y/Z^A ln ^ ont . Eq. (B7l is used in section [v| while calculating the noise back-action factor 

n. 

The validity of this crossover to quantum description lies in the fact that in case of parametric interaction, the 
difference between the classical and quantum evolution vanishes when the number of photons in the line is large or the 
coupling of the system to reservoir is wealP^. We can then regard the quantum fluctuations to be driven by classical 
random fields, obeying classical equations of motion. 
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APPENDIX C 



Table of Symbols 


A 


Reduced detuning of signal from u g {= " s r " 8 j 


Ej 


Josephson energy 


F 


Effective pump strength 


F B a 


Back-action corrected pump strength 


G 


Power gain of the amplifier 


r 


Damping rate {2W0) 


7 


Effective nonlinearity parameter ( = r , , m \\ , , I 


©eff 


hnective temperature at the input port (= ^y) 


V 


Back-action index (= 7 • B c ff) 


In 


Noise current 


Irf 


RF drive current 


Lj 


Josephson inductance 


$ 


Flux across the Josephson Junction 


$0 


Mux quantum (= 




Dimcnsionlcss flux variable 


n 


li UloL UCJiVjiS. Civ, LilUll 1 ClG L \JL 




T-? penn s nt frpni ion pit - at TAQonriQAn nQPillsiiTn'" 




Piittiti TrpniiPTir'v 1 




Pump frequency 2 




Ghost freaucnev (— Ul + U2 ) 
irequency ^ 2 ; 


o 


neciuceo cierunmg 01 ujg nom ujq i — p i 


r 


Reflection coefficient 


s 


Interconversion gain coefficient 




Photon number spectral density 


T 


Black body temperature at the input (signal and idler ports) 


T N 


Noise temperature of the amplifier 


X[w] 


Amplitude of noise component at generic frequency u> 




Steady state amplitude phasor 


Y 


Amplitude phasor for pump 1 


Z 


Amplitude phasor for pump 2 
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